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Abstract 

Superconducting mesoscopic devices in magnetic fields present novel prop- 
erties which can only be accounted for by both the quantum confinement of 
the Cooper pairs and by the interaction between the magnetic-field-induced 
vortices. Sub-micrometer disks, much the same as their semiconductor coun- 
terparts known as quantum dots, are being subject to experimental inves- 
tigation by measuring their conducting properties and, more recently, their 
magnetization by using state-of-the-art ballistic Hall magnetometry. In this 
work I review the main results obtained in these two types of experiments 
as well as the current theoretical developments which are contributing to our 
understanding of the superconducting condensate in these systems. 
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I. INTRODUCTION 



The physics of mesoscopic systems has reached its maturity within the semiconductor 
community over the past few years Q. A similar revolution is nowadays taking place in 
superconductivity where superconducting device miniaturization is moving forward at a fast 
pace ||. In both communities the concept of mesoscopics is directly linked to that of low- 
dimensionality. A system is called low- dimensional when at least one of the dimensions in 
which its carriers live becomes comparable to some relevant length scale intrinsically asso- 
ciated to bulk properties of those carriers. In semiconductors the most relevant length scale 
one can think of is the Fermi wavelength Af- In superconductors there are two fundamental 
length scales. One is the penetration length A over which an externally applied magnetic 
field H is screened into the superconducting condensate. The other one is the coherence 
length £ which, in a simple-minded picture, can be thought of as the size of the Cooper 
pairs and which limits the distance over which the superconducting order parameter can 
vary appreciably. When one of the dimensions of a semiconductor becomes comparable to 
Af, the other two being much larger, a two-dimensional electron gas (2DEG) is formed [|TJ. 
Similarly, typical superconducting thin films can have a width comparable to either A 
or £, or to both of them. If, in semiconductors, two out of the three dimensions become 
comparable to Ap one forms a quantum wire |T|. A thin strip || can be thought of as the 
superconducting counterpart. When the three dimensions become comparable to either Af 
in semiconductors and to A and/or £ in supercondutors one speaks of 'zero- dimensional' 
systems called quantum dots (§-|| and thin disks respectively. When the dimensions 

of the superconducting device become much smaller than £ (as can be the case in very small 
metallic grains) the very notion of superconductivity needs to be revised |Uj. The above 
classification, far from being unique and rigorous, only pretends to serve as a guide for the 
non-specialized reader. 

It has been known for a long time that, in addition to the well- understood type-I su- 
perconductors where k = A/£ < l/y/2, superconductors with k > l/v2 are also possible 
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(high-T c materials are the most exotic example of them ||). This seemingly inocent relation 
between the two fundamental lengths gives rise to type-II superconductivity which is charac- 
terized by the appearance of vortices in the superconducting condensate at values of H lying 
between a lower critical field H c \ and an upper critical field H C 2. An isolated vortex extends 
over a distance of the order of A and exhibits a normal core of radius of the order of £ 0. A 
macroscopic number of vortices generically form a triangular lattice |l]J which presents two 
well- characterized regimes or states: (i) a dilute vortex state (dVS) at low fields in which A 
governs the inter-vortex distance and (ii) a dense vortex state (DVS) at high fields in which 
the fundamental length scale is given by £ since the vortices are closely packed and their 
cores overlap strongly. In extreme type-II superconductors, characterized by k ^> 1, the 
crossover from the dVS to the DVS takes place at rather low fields H pa 0.3H c2 |T3| which 
gives an idea of the relative importance of both regimes. 

The dVS has been thoroughly studied over the years in all possible low-dimensional sys- 
tems, both theoretical and experimentally. Many experimental techniques are sensitive to 
the vortex lattice in this regime, mainly due to the strong spatial modulation of the asso- 
ciated magnetic induction ]r2| . On the theory part, vortices in the dVS can be treated as 
classical, logarithmically-interacting, point-like particles which renders the calculation of the 
vortex structure free energy rather feasible in many different geometries [|13H T5|. By con- 
trast, the structure of the DVS in nanostructures is difficult to unveil, partially due to the 
'invisibility' of this state to the usual experimental techniques and to the difficulty in mini- 
mizing the free energy associated to this highly compacted vortex state when interfaces are 
present. Motivated by very recent advances in transport [0] and magnetometry techniques 
||, which we review in Sec. II, we focus in this work on the theoretical study of the struc- 
tural and magnetic properties of the DVS in superconducting thin disks or quantum dots. 
In Sec. Ill we show how to do this within the traditional framework of the Ginzburg-Landau 
functional in a rather analytical fashion. Finally, in Sec. IV we present our conclusions. 
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II. TWO RELEVANT EXPERIMENTS IN SUPERCONDUCTING MESOSCOPIC 

DISKS 



One the most remarkable achievements in low-dimensional semiconductor physics has 
been the fabrication of the quantum dot or single-electron transistor where the electrons 
travel through the system one at a time and where the number of them present in the device 
can be tuned at will even down to a single electron [[|||. These systems have been also 
given the name of artificial atoms since their generic properties are determined by both the 
quantum confinement and the interaction among the electrons much the same as in real 
atoms. Many analogies can be drawn between these artificial atoms and superconducting 
mesoscopic disks in perpendicular magnetic fields. Cooper pairs in disks can also experience 
the effects of quantum confinement and this becomes visible in the superconducting-normal 
phase boundary 0. As far as the interaction is concerned the role of the electron is now 
played to some extent, not by the Cooper pair, but by another fundamental entity: The 
vortex. When the dimensions of the disk are comparable to A in the dVS or to £ in the DVS 
only few vortices can coexist in the disk. In contrast to the usual triangular arrangement in 
bulk, complex and unique vortex structures are expected to occur due to the competition 
between the geometrical confinement and the vortex-vortex interaction. 

Transport experiments contributed in a decisive way to unveil the electronic structure 
of artificial atoms ||. Similarly, conductance measurements [|7| in individual mesoscopic 
disks gave us the first experimental evidence of the structure of the order parameter in these 
systems. Based on the onset of the disk resistance, oscillations of the critical temperature 
T c as a function of the external magnetic field were measured and correctly accounted for 
by the quantization of the angular momentum L of the Cooper pair wavefunction. In more 
traditional words, they were observing transitions between giant vortex states with a different 
number L of fluxoids. These oscillations were clearly reminiscent of those observed in the 
famous Little-Parks experiment |16| . 



Transport measurements are necessarily bound to give only information about the crit- 
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ical phase boundaries since they are based on measuring the resistance of the non-fully 
superconducting state of the disk. Hall magnetometry ||, on the other hand, is revealing 
itself as a powerful tool for obtaining information of the order parameter away from the 
supercondutor-normal phase boundary Thin superconducting Al disks are placed on top of 
a series of typical Hall probes created in a 2DEG. In the ballistic regime the Hall resistance is 
directly proportional to the average value of the magnetic field through the junction which, 
in turn, is determined by the magnetic state of the disk lying right above the junction. 
Magnetization (M) curves for different disk sizes can thus be obtained to exhibit a variety 
of unexpected phase transitions. For very small disks a first-order transition at a critical 
H c , which is what one expects for a type-I superconductor as Al, is absent. For larger sizes 
this transition appears, but, on increasing further the disk radius the familiar second-order 
phase transition for a type-II superconductor, namely, an steady decrease of the magnetiza- 
tion beyond a certain critical H&, becomes notorious. On top of this steady decrease, the 
magnetization exhibits many jumps corresponding to first-order transitions which present 
an irregular and decreasing height as a function of H. All this seems to indicate that the 
quantum confinement is also patent in this experiment and that the Al disks are behaving 
like type-II superconductors rather than type-I, possibly due to the expected enhancement 
of the effective magnetic penetration length in such a geometry. 

In the next section we show how many of these unexpected features can be obtained and 
explained within the traditional, phenomenological Ginzburg-Landau theory. 

III. GINZBURG-LANDAU FUNCTIONAL APPROACH TO THE DVS IN 
SUPERCONDUCTING MESOSCOPIC DISKS 

The theoretical analysis of the peculiar supercondutor-normal phase boundary measured 
in Ref. |7[] does not present any difficulty since it simply implies solving the linearized 
Ginzburg-Landau equations [|?yiT|-pOH . The theoretical efforts to calculate the properties of 
the superconducting condensate in a disk well into the superconducting phase have been 
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mostly restricted to the dVS []14| , |15[ . As far as the DVS is concerned, recent work is being 
done in the direction of solving numerically the Ginzburg-Landau equations either under 
the simplifying assumption of an order parameter with a well-defined L [pO| , pT| or without 



any symmetry restrictions P J2^J23|1 . For type-II disks, the assumption of an order parameter 
with axial symmetry (well-defined L) does not hold. More precisely, it is only expected to 
hold in the Meissner state, which is associated to an L = order parameter, and above H C 2 



24j where, close to the surface of the disk, the superconductivity can survive up to a higher 



critical field H c3 [|7l,|i7HH|l . This surface superconductivity is referred to as a condensed 



state of vortices or giant vortex [20 — 23 1 . For H c \ < H < H C 2 it was argued pC§ and recently 



shown in numerical simulations |8|j23|| that the order parameter can form complex structures 



of single-fluxoid vortices, i.e., a "budding" Abrikosov lattice. From our more analytic stand- 
point the appearance of these structures necessarily implies an order parameter without a 
well-defined angular momentum L. Recent work in which this condition is explicitly taken 



into account has been presented in Ref. and also in Ref. . Part of the results presented 
here have already appeared in Ref. [[J. 

We start from the traditional Ginzburg-Landau functional for the Gibbs free energy of 
the superconducting state 



G s = G a + J dr 



a|*(r)| 2 + f|*(r)| 4 + 
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where G n is the Gibbs free energy of the normal state and [— iTiS7 — e*A(r)/c] 2 /2m* is the 
kinetic energy operator for Cooper pairs of charge e* = 2e and mass m* = 2m in a vector 
potential A(r) which is associated with the magnetic induction h(r). The parameters a 
and (3 have the usual meaning 0. Before proceeding any further we must stress a not fully 
appreciated fact: Even for small values of k (~ 1), the magnetic induction is weakly varying 
in space down to fairly low fields (H w 0.5H C 2) |I2| . Thus, it is a very good approximation 
to consider a uniform magnetic induction [h(r) = B] down to H w 0.5H C 2 and to expand 
the order parameter in the lowest Landau level (LLL): 
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*(r) = £c L — =e-* Lfl $ L (r). (2) 
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In this expansion = \Ci\e l ^ L are complex coefficients and jj^ e ~ %Ld ^ are normalized 
eigenfunctions of the kinetic energy operator in Eq. ([!]) plus the boundary conditions of zero 
current through the surface of the disk . (We only consider disk thicknesses smaller than 



the coherence length so that the order parameter can be taken constant in the direction 
of the field.) Strictly speaking, these eigenfunctions are not LLL eigenstates. However, 
their radial part $L(r), which may be found numerically, is nodeless and coincides with 
the radial function of the symmetric LLL eigenfunctions for small L. Figure [l] shows the 
Cooper pair band structure, i.e., the corresponding eigenvalues at different values of the 
external field for a disk of radius R = 8£ and k = oo (this implies H = B). The horizontal 
lines represent —a which can be thought of as the chemical potential for the Cooper pairs. 
The formation of the dispersionless LLL in the center of the disk as the magnetic length 



'e*Ti/cB becomes smaller than R is clearly visible. The bending of the band close to 
the surface is a consequence of the boundary condition and allows the nucleation of surface 
superconductivity (or the formation of a giant vortex) even when the bulk remains in the 
normal state (H > H C 2). 

The expansion (Q) captures both the simplicity of the macrovortex (above H C 2) when only 
one Cl is expected to be different from zero and the full complexity of the order parameter 
(below H C 2) when several components or harmonics must participate. Direct substitution 
of such an expansion into Eq. [I] and subsequent numerical minimization of the resulting 
expression is a daunting task bound to fail due to the large number of unknown variables 
involved. Instead, it is preferable to consider expansions in restricted sets {Li, L 2 , . . . , L^} 
of few iV components where L\ < L2 < . . . < Ljy. The difference between the Gibbs free 
energies of the normal and superconducting phases takes the following form for each set: 

N 1 

G s - G n = 5>[1 - Be Li (B)]\C Li \ 2 + -a 2 n 2 BR 2 x 
i=i 4 

' N N 
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where G s — G n and a are given in units of H 2 2 V/8n (V is the volume of the disk), €l(B) 
is given in units of the lowest Landau level energy Tiu c /2 (with oo c = e*B/m*c), R is in 
units of £, and B and H are given in units of H c2 . The terms proportional to a contain 
the condensation and kinetic energy of the Cooper pairs. All the other terms, which are 
proportional to a 2 , account for the "interaction" between Cooper pairs. There appear four 
types of these terms: (i) those proportional to Il{B) = / dr r reflecting the interac- 
tion between Cooper pairs occupying the same quantum state L, (ii) those proportional to 
iLiLj(B) = / dr r reflecting the interaction between Cooper pairs occupying differ- 

ent quantum states, and (iii) the ones proportional to l^ 3 ^ k [B) = / dr r $^$1 $£, fc and 
proportional to Il.l^lXB) = J dr r $ Li $ Lj $ Lfc $ L; which, along with the phases <f> L , are 
responsible for the spatial correlation between vortices and the detailed structure of the DVS 
(see below). The non-linear dependence on B of these integrals [as well as that of €l(B)] 
comes from the existence of the disk surface. 

In order to find the minimum Gibbs free energy for a given set we have to minimize 
with respect to the moduli \Cl x |, . . . , IC^I, the phases (pn, ■ ■ ■ , 4>l n of the coefficients, and 
with respect to B. The minimum-energy set of components is picked up at the end. The 
advantage of doing this selective minimization resides in our expectation that a small number 
of components will suffice to describe the order parameter for disks with radii of few coherent 
lengths. As illustrative and relevant examples we consider in the detail the solutions with 
one and two components. For N — 1 the energy functional is invariant with respect to the 
phase of the only coefficient so one can minimize analytically with respect to \Cl\ 2 to obtain 



where we have dropped the implicit 5-dependences. Finally, the minimal value of B and 
the minimum Gibbs free energy for each L must be found numerically. The 2-component 
solutions {Li,L 2 } can be dealt with in a similar way. The energy functional is invariant 
with respect to the phases so one can minimize analytically with respect to |CxJ 2 and |Cx 2 | 2 
to obtain 

r r [I- Be Ll fl L2 + [l- Bt L2 fl Ll -A\l- Be Ll ]\l- Be L2 ]I LlL2 . 

Gr s -G n = ^2um\T t 17 1 V[B-H) (5 J 

and numerically with respect to B to obtain the minimum energy for the given pair of 
components. The superconducting condensate looks generically like an (L 2 — Li)-vortex ring. 
It is important to notice that a disappears from the final expressions (f|) and @ which leaves 
us with k as the only adjustable parameter when comparing with experiments. (This is also 
true for the more complex cases discussed below). For 3-component solutions (two vortex 
rings) the energy functional is still invariant with respect to all the three phases whenever 
L 1 + L 3 7^ 2L 2 , and, once again, the minimization with respect to (C^J 2 , |Cl 2 | 2 , and |Cx 3 | 2 
can be done analytically. However, if L\ + L3 = 2L 2 , the two rings have the same number of 
vortices and their relative angular positions come into play through the term depending on 
the phases. There is, however, an obvious choice for these phases: 0^ =0, <pL 2 — 0, 0l 3 = 
This choice gives a negative contribution to the free energy which reflects a lock-in position 
between the vortex rings. In general, the minimum-energy solutions for disks are expected to 
have strongly overlapping components which invalidates any perturbative treatment of the 
terms that contain the phases ||. Moreover, unlike simpler geometries ||, there is no direct 
connection between number of components in which we expand the order parameter and 
number of vortices. This prompts us to seek solutions through numerical minimization with 
respect to the moduli and B for the 3-component cases just mentioned, and, for N > 3, with 
respect to the moduli, the phases, and B whenever the terms involving phases are present. 
For the disk sizes like the ones used in the experiment of Ref. || 2 and 3-component solutions 
suffice to capture the relevant physics. 

Figure ||] shows the magnetization as a function of H for two disks with k = 2. The com- 
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ponents or harmonics with significant weight corresponding to the minimum-energy order 
parameter at each magnetization step are also shown. For the smaller disk of radius R = 4£ 
we obtain a series of first-order transitions. Here all the minimal solutions correspond to 
giant vortices which contain L fluxoids. Whenever L changes by one the magnetization 
presents a (non-quantized) jump whose magnitude evolves monotonously with L. For the 
larger R = 7£ disk N > 1 solutions appear below H C 2. In this case, allowing more compo- 
nents in the expansion of the order parameter has a fundamental effect: It splits the giant 
vortex into a complex vortex glass structure of many single-fluxoid vortices [see Figs. ||(a) 
to (c)]. This glass structure reflects in the magnetization curves by changing the regular 
evolution of the magnitude of the jumps into an irregular one. Whenever a vortex is added 
or removed from the disk the symmetry of the new vortex configuration changes which, in 
turn, expels the field in a different way. The total number of vortices is always given by 
the largest which does not depend on the value of L of the other harmonics: Only the 
internal arrangement of vortices does. There usually exist configuration switches for a given 
number of vortices (i.e., for a given L/v)> but these changes do not reflect in the magnetiza- 
tion (see Fig |2|). On top of the first-order transitions the overall slope in the magnetization 
clearly changes at H c2 , i.e., at the transition between the giant vortex and the vortex glass 
structures. This transition is reminiscent of the second-order transition at H C 2 for bulk 
samples where M vanishes. 

Finally, we point out that the magnetization measured by Geim et al. in Ref. || presents 
features that are in good agreement with our results despite of the fact that the disks are 
made out of a strong type-I material as Al. It is well known that the bulk value of k can 



increase in the plate geometry pqj . To compare with the experiment, we simulate this 
fact by using a higher value of k than the nominal one. In Fig. ^ we show the data for a 
disk of nominal radius R = 5£ and thickness d = 0.6£. We have obtained a reasonable good 
agreement in the number and magnitude of the jumps, and overall shape of the magnetization 
using R « 5£ and k « 1 (the dotted line is a good example). This is consistent with an 
effective penetration length longer than expected and, possibly, with a coherence length 
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shorter than the bulk nominal one. Although, due to the smallness of the disk, it is difficult to 
point at a second-order phase transition, the non-monotonous evolution of the magnitude of 
the magnetization jumps is notorious over a large range of fields which, as we have shown, is 
an indication of the formation of vortex glass structures. In our approximation the magnetic 
induction is uniform in space, but, given the good agreement with the experimental curve, 
this does not seem to be an important restriction. 

IV. CONCLUSIONS 

After reviewing some of the state-of-the-art experiments in mesoscopic superconductivity 
we have shown how to calculate the dense vortex matter structure and associated magneti- 
zation for type-II superconducting mesoscopic disks. We have found that the magnetization 
exhibits generically a first-order phase transition whenever the number of vortices changes by 
one with H. It also presents two well-defined regimes: A non-monotonous evolution of the 
magnitude of the jumps signals the presence of a vortex glass structure which is separated 
by a second-order phase transition at H C 2 from a condensed state of vortices (giant vortex) 
where the magnitude of the jumps changes monotonously. We have compared our results 
with the Hall magnetometry measurements in Ref. [g] and claimed that the magnetization 
exhibits clear traces of the presence of these vortex glass states. 
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FIG. 1. Band structure and chemical potential (—a) of the Cooper pairs at different values of 
H (in units of H C 2) for a disk of radius R = 8£ and k = oo. Notice that, for fields above i? C 2, the 
lowest Landau level can lie below the chemical potential only close to the surface of the disk. This 
gives rise to the formation of a giant vortex with an L around the bottom of the band. 
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FIG. 2. Magnetization as a function of H for two disks of k = 2 and radius (a) R = 4£ and (b) 
R = 7£. The sets of harmonics with significant weight corresponding to the expansion of the order 
parameter at each magnetization step are also shown. 




FIG. 3. Superconducting condensate in the DVS for an R = 7£, k = 2 disk at (a) H = 0.5H C 2, 
(b) H = 0.7H C 2, (c) H = 0.9i7 C 2, and (d) H = 1.1H C 2. Below H C 2 the density associated to the 
order parameter presents a vortex glass structure which disappears above H C 2 into a condensate of 
vortices or giant vortex. Notice the appearance in (b) of two-fluxoid vortices even below H C 2. 
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FIG. 4. Comparison between experimental data (reproduced from Refs. [21]) and theory using 
R = 5.25£ and k = 1.2. Similar considerations as in Refs. [21] have been followed for the adjustment 
of the theoretical curve. 
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